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I summarize the current state of solar neutrino research and then give my answer to the question: What should 
we do next? 



1. Introduction 

The reader who is familiar with solar neutrino 
research may wish to skip directly to the last sec- 
tion entitled: What Next? 

Solar neutrinos have been detected experimen- 
tally with fluxes and energies that are quali- 
tatively consistent with solar models that are 
constructed assuming that the sun shines by 
nuclear fusion reactions. The first experimen- 
tal result, obtained by Ray Davis and his col- 
laborators in 1968, has now been con- 
firmed by four other beautiful experiments, 
Kamiokande, @ SAGE, g GALLEX, @ and 
SuperKamiokande. j(| The observation of solar 
neutrinos with approximately the predicted en- 
ergies and fluxes establishes empirically the the- 
ory |?J that main sequence stars derive their en- 
ergy from nuclear fusion reactions in their inte- 
riors and has inaugurated what we all hope will 
be a flourishing field of observational neutrino as- 
tronomy. 

Although the calculated neutrino fluxes depend 
upon high powers of the central temperature of 
the solar model, the experiments and the solar 
model theory are so precise that persistent quan- 
titative discrepancies have existed between the 
model predictions and the solar model calcula- 
tions for over thirty years §-|h]] 

Important experiments are underway that 
will provide diagnostic information about the 
physical properties of neutrinos that are cre- 
ated in the center of the sun and detected 
on earth in really long baseline experiments. 
At this workshop, we will hear discussions 
of the SuperKamiokande, SNO, BOREXINO, 



HELLAZ, HERON, ICARUS, LENS, and Kam- 
LAND experiments. 

I will discuss predictions of the combined stan- 
dard model in the main part of this review. By 
'combined' standard model, I mean the predic- 
tions of the standard solar model and the pre- 
dictions of the standard electroweak model. We 
need a solar model to tell us how many neutri- 
nos of what energy are produced in the sun and 
we need an electroweak theory to tell us how the 
number and flavor content of the neutrinos are 
changed as they make their way from the center 
of the sun to detectors on earth. For all practical 
purposes, the standard electroweak model states 
that nothing happens to solar neutrinos after they 
are created in the deep interior of the sun. Us- 
ing standard electroweak theory and fluxes from 
the standard solar model, one can calculate the 
rates of neutrino interactions in different terres- 
trial detectors with a variety of energy sensitivi- 
ties. The combined standard model also predicts 
that the energy spectrum from a given neutrino 
source should be the same for neutrinos produced 
in terrestrial laboratories and in the sun and that 
there should not be measurable time-dependences 
(other than the seasonal dependence caused by 
the earth's orbit around the sun). The spectral 
and temporal departures from standard model ex- 
pectations are expected to be small in all cur- 
rently operating experiments |ll}| and have not 
yet yielded definitive results. Therefore, I will 
concentrate here on inferences that can be drawn 
by comparing the total rates observed in solar 
neutrino experiments with the combined stan- 
dard model predictions. 

I will begin by reviewing in Section || the quan- 



2 



titative predictions of the combined standard so- 
lar model and then describe in Section || the three 
solar neutrino problems that are established by 
the chlorine, Kamiokande, SAGE, GALLEX, and 
SuperKamiokande experiments. In Section ||, I 
detail the uncertainties in the standard model 
predictions and then show in Section || that he- 
lioseismological measurements indicate that the 
standard solar model predictions are accurate for 
our purposes. In Section ||, I discuss the im- 
plications for solar neutrino research of the pre- 
cise agreement between helioseismological mea- 
surements and the predictions of standard solar 
models. Next, ignoring all knowledge of the sun, I 
cite analyses in Section ^| that show that one can- 
not fit the existing experimental data with neu- 
trino fluxes that are arbitrary parameters, unless 
one invokes new physics to change the shape or 
flavor content of the neutrino energy spectrum. 
I summarize in Section ^ the characteristics of 
the best-fitting neutrino oscillation descriptions 
of the experimental data. Finally, I will discuss 
and summarize the results in Section |§|. 

If you want to obtain numerical data or subrou- 
tines that are discussed in this talk, or to see rele- 
vant background information, you can copy them 
from my Web site: bttp : / / www. sns . ias . cdu/ ~ j nt . 

Before we begin the detailed discussion, I want 
to make just a brief historical diversion. Nearly 
all of the current interest in solar neutrinos cen- 
ters around the opportunity to use the sun as a 
neutrino source in a very long baseline oscillation 
experiment. In preparing for a talk in honor of 
Fred Reines a few months ago, I ran across a long 
forgotten 1972 letter from Bruno Pontecorvo, the 
originator of the hypothesis that oscillations may 
be observed in solar neutrino experiments. For 
your interest, I enclose a reproduction of this let- 
ter in Fig. |l]. 

For the benefit of neutrino pioneers of today, it 
is perhaps worth remarking that Ray Davis and I 
never considered the possibility that solar neutri- 
nos could be used to learn more about neutrinos 
when, in the early 1960's, we were first analyzing 
the potentialities of a practical chlorine experi- 
ment. We sold the experiment as a fundamental 
test of the hypothesis that the sun shines by nu- 
clear fusion reactions in its interior. Only after 
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Figure 1. Letter from Bruno Pontecorvo in 1972. 



the first results of the chlorine experiment showed 
in 1968 a conflict with the solar model calcula- 
tions and Gribov and Pontecorvo published their 
epochal 1969 paper on vacuum oscillations of so- 
lar neutrinos did we begin to consider the possi- 
bility that solar neutrinos might tell us something 
new about particle physics. Maybe there are pre- 
viously unimagined physics treasures to be dis- 
covered in future neutrino experiments. 

2. Standard Model Predictions 

Table @ gives the neutrino fluxes and their 
uncertainties for our best standard solar model, 
hereafter BP98. |1(J Figure || shows the pre- 
dicted neutrino fluxes from the dominant p-p fu- 
sion chain. 

The BP98 solar model includes diffusion of 
heavy elements and helium, makes use of the nu- 
clear reaction rates recommended by the expert 
workshop held at the Institute of Nuclear The- 
ory, recent (1996) Livermore OPAL radiative 
opacities, jl3| the OPAL equation of state, 
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Table 1 

Standard Model Predictions (BP98): solar neu- 
trino fluxes and neutrino capture rates, with lc 
uncertainties from all sources (combined quadrat - 
ically). 



Gallium [■ 



Chlorine i 



SuperK 



Source 



Flux 
'to 10 cm^s- 1 ) 



CI Ga 
(SNU) (SNU) 



PP 

pep 

hep 

7 Be 

8 B 

13 N 
15 o 

17 F 
Total 



5.94 (l.OOtHx) 



i-00±8;^) 



1.39 x 10~ 2 
2.10 x 1(T 7 
4.80 x 10- 1 (l.OOtoil) 
5.15 x 10- 4 (l.00t°;^) 

6.05 x io- 2 (l.ootgil) 

5.32 x 10~ 2 (l.OOt^; 22 ) 

6.33 x 10~ 4 (l.OOt^ 2 ) 



0.0 

0.2 

0.0 

1.15 

5.9 

0.1 

0.4 

0.0 



69.6 

2.8 

0.0 

34.4 

12.4 

3.7 

6.0 

0.1 



7.7±i;o I29lj? 



and electron and ion screening as determined by 
the recent density matrix calculation. p5| , ^6[ The 
neutrino absorption cross sections that are used 
in constructing Table [l] are the most accurate val- 
ues available [ p"7|]l8| and include, where appropri- 
ate, the thermal energy of fusing solar ions and 
improved nuclear and atomic data. The valid- 
ity of the absorption cross sections has recently 
been confirmed experimentally using intense ra- 
dioactive sources of 51 Cr. The ratio, R, of the 
capture rate measured (in GALLEX and SAGE) 
to the calculated 51 Cr capture rate is R = 0.95 ± 
0.07 (exp) + to°oi (theory) and was discussed 
extensively at Neutrino 98 by Gavrin and by 
Kirsten. The neutrino-electron scattering cross 
sections, used in interpreting the Kamiokande 
and SuperKamiokande experiments, now include 
electroweak radiative corrections. []l9f 

Figure shows for the chlorine experiment all 
the predicted rates and the estimated uncertain- 
ties (lcr) published by my colleagues and myself 
since the first measurement by Ray Davis and 
his colleagues in 1968. This figure should give 
you some feeling for the robustness of the solar 
model calculations. Many hundreds and prob- 
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Figure 2. The energy Spectrum of neutrinos from 
the pp chain of interactions in the Sun, as pre- 
dicted by the standard solar model. Neutrino 
fluxes from continuum sources (such as pp and 
8B) are given in the units of counts per cm2 per 
second. The pp chain is responsible for more 
than 98% of the energy generation in the stan- 
dard solar model. Neutrinos produced in the 
carbon-nitrogen-oxygen CNO chain are not im- 
portant energetically and are difficult to detect 
experimentally. The arrows at the top of the fig- 
ure indicate the energy thresholds for the ongoing 
neutrino experiments. 



ably thousands of researchers have, over three 
decades, made great improvements in the input 
data for the solar models, including nuclear cross 
sections, neutrino cross sections, measured ele- 
ment abundances on the surface of the sun, the 
solar luminosity, the stellar radiative opacity, and 
the stellar equation of state. Nevertheless, the 
most accurate predictions of today are essentially 
the same as they were in 1968 (although now they 
can be made with much greater confidence). For 
the gallium experiments, the neutrino fluxes pre- 
dicted by standard solar models, corrected for dif- 
fusion, have been in the range 120 SNU to 141 
SNU since 1968. (l7) A SNU is a convenient unit 
with which to describe the measured rates of so- 
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Predicted 37 C1 Rate vs. Time 
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Figure 3. The predictions of John Bah- 
call and his collaborators of neutrino capture 
rates in the 37 C1 experiment are shown as 
a function of the date of publication (since 
the first experimental report in 1968. Q) The 
event rate SNU is a convenient product of 
neutrino flux times interaction cross section, 
10~ 36 interactions per target atom per sec. The 
format is from Figure 1.2 of the book Neutrino 
Astrophysics. |j] The predictions have been up- 
dated through 1998. 



lar neutrino experiments: 10~ 36 interactions per 
target atom per second. 

There are three reasons that the theoretical cal- 
culations of neutrino fluxes are robust: 1) the 
availability of precision measurements and preci- 
sion calculations of input data; 2) the connection 
between neutrino fluxes and the measured solar 
luminosity; and 3) the measurement of the he- 
lioseismological frequencies of the solar pressure- 
mode (p-mode) eigenfrequencies. I have discussed 
these reasons in detail in another talk. J|(| 

Figure ^ displays the calculated 7 Be and 8 B 
neutrino fluxes for all 19 standard solar models 
which have been published in the last 10 years 
in refereed science journals. The fluxes are nor- 
malized by dividing each published value by the 
flux from the BP98 solar model; [[l0| the abscissa 
is the normalized 8 B flux and the ordinate is the 
normalized 7 Be neutrino flux. The rectangular 



box shows the estimated 3er uncertainties in the 
predictions of the BP98 solar model. 

All of the solar model results from different 
groups fall within the estimated 3<r uncertain- 
ties in the BP98 analysis (with the exception 
of the Dar-Shaviv model whose results have not 
been reproduced by other groups). This agree- 
ment demonstrates the robustness of the predic- 
tions since the calculations use different computer 
codes (which achieve varying degrees of precision) 
and involve a variety of choices for the nuclear 
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Figure 4. Predictions of standard solar mod- 
els since 1988. This figure, which is Fig. 1 of 
Bahcall, Krastev and Smirnov (1998) JO]], shows 
the predictions of 19 standard solar models in 
the plane defined by the 7 Be and 8 B neutrino 
fluxes. The abbreviations that are used in the fig- 
ure to identify different solar models are defined 
in the bibliographical item, Ref. |2l|]. The figure 
includes all standard solar models with which I 
am familiar that were published in refereed jour- 
nals in the decade 1988-1998. All of the fluxes 
are normalized to the predictions of the Bahcall- 
Pinsonneault 1998 solar model, BP98. fo| The 
rectangular error box defines the 3cr error range 
of the BP98 fluxes. The best-fit 7 Be neutrino flux 
is negative. At the 99% C.L., there is no solu- 
tion jll|] with all positive neutrino fluxes (see dis- 
cussion in Section ^) . All of the standard model 
solutions lie far from the best-fit solution, even 
far from the 3(7 contour. 



5 



parameters, the equation of state, the stellar ra- 
diative opacity, the initial heavy element abun- 
dances, and the physical processes that are in- 
cluded. 

The largest contributions to the dispersion in 
values in Figure f| are due to the choice of the nor- 
malization for Sn (the production cross-section 
factor for 8 B neutrinos) and the inclusion, or non- 
inclusion, of clement diffusion in the stellar evolu- 
tion codes. The effect in the plane of Fig. || of the 
normalization of Sn is shown by the difference 
between the point for BP98 (f.O,f.O), which was 
computed using the most recent recommended 
normalization, Jl^j and the point at (f.f8,f.O) 
which corresponds to the BP98 result with the 
earlier (CalTech) normalization, p^ ] 

Helioseismological-observations 
have shown 



Total Rates: Standard Model vs. Experiment 

Bahcall-Pinsonneault 98 



that element diffusion is oc- 
curring and must be included in solar models, so 
that the most recent models shown in Fig. [| now 
all include helium and heavy element diffusion. 
By comparing a large number of earlier models, 
it was shown that all published standard solar 
models give the same results for solar neutrino 
fluxes to an accuracy of better than 10% if the 
same input parameters and physical processes are 
included. @,|| 

Bahcall, Krastev, and Smirnov |[l] have com- 
pared the observed rates with the calculated, 
standard model values, combining quadratically 
the theoretical solar model and experimental un- 
certainties, as well as the uncertainties in the 
neutrino cross sections. Since the GALLEX and 
SAGE experiments measure the same quantity, 
we treat the weighted average rate in gallium as 
one experimental number. We adopt the Super- 
Kamiokande measurement as the most precise di- 
rect determination of the higher-energy 8 B neu- 
trino flux. 

Using the predicted fluxes from the BP98 
model, the % 2 for the fit to the three experimental 
rates (chlorine, gallium, and SuperKamiokandc, 
see Fig. ||) is 



Xg SM (3 experimental rates) = 61 



(1) 



The result given in Eq. (0), which is approxi- 
mately equivalent to a 20a discrepancy, is a quan- 
titative expression of the fact that the standard 
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Figure 5. Comparison of measured rates and 
standard-model predictions for five solar neutrino 
experiments. The unit for the radiochemical 
experiments (chlorine and gallium) is SNU (see 
Fig. U for a definition); the unit for the water- 
Cerenkov experiments (Kamiokande and Super- 
Kamiokande) is the rate predicted by the stan- 
dard solar model plus standard electroweak the- 
ory. M 



model predictions do not fit the observed solar 
neutrino measurements. 

3. Three Solar Neutrino Problems 

I will now compare the predictions of the com- 
bined standard model with the results of the op- 
erating solar neutrino experiments. 

We will see that this comparison leads to three 
different discrepancies between the calculations 
and the observations, which I will refer to as the 
three solar neutrino problems. 

Figure || shows the measured and the calculated 
event rates in the five ongoing solar neutrino ex- 
periments. This figure reveals three discrepancies 
between the experimental results and the expec- 
tations based upon the combined standard model. 
As we shall see, only the first of these discrepan- 
cies depends in an important way upon the pre- 
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dictions of the standard solar model. 

3.1. Calculated versus Observed Absolute 
Rate 

The first solar neutrino experiment to be per- 
formed was the chlorine radiochemical experi- 
ment, H which detects electron-type neutrinos 
that are more energetic than 0.81 MeV. After 
more than a quarter of a century of operation 
of this experiment, the measured event rate is 
2.56 ± 0.23 SNU, which is a factor of three less 
than is predicted by the most detailed theoret- 
ical calculations, 7.7t\l SNU. Jl§ Most of the 
predicted rate in the chlorine experiment is from 
the rare, high-energy 8 B neutrinos, although the 
7 Be neutrinos are also expected to contribute sig- 
nificantly. According to standard model calcula- 
tions, the pep neutrinos and the CNO neutrinos 
(for simplicity not discussed here) are expected 
to contribute less than 1 SNU to the total event 
rate. 

This discrepancy between the calculations and 
the observations for the chlorine experiment was, 
for more than two decades, the only solar neutrino 
problem. I shall refer to the chlorine disagreement 
as the "first" solar neutrino problem. 

3.2. Incompatibility of Chlorine and Water 
Experiments 

The second solar neutrino problem results from 
a comparison of the measured event rates in the 
chlorine experiment and in the Japanese pure- 
water experiments, Kamiokande Q and Super- 
Kamiokande. || The water experiments detect 
higher-energy neutrinos, most easily above 7 
MeV, by by observing the Cerenkov radiation 
from neutrino-electron scattering: v + e — > 
v' + e'. According to the standard solar model, 8 B 
beta decay, and possibly the hep reaction, ||(J are 
the only important source of these higher-energy 
neutrinos. 

The Kamiokande and SuperKamiokande exper- 
iments show that the observed neutrinos come 
from the sun. The electrons that are scattered 
by the incoming neutrinos recoil predominantly 
in the direction of the sun-earth vector; the rel- 
ativistic electrons are observed by the Cerenkov 
radiation they produce in the water detector. In 



addition, the water Cerenkov experiments mea- 
sure the energies of individual scattered electrons 
and therefore provide information about the en- 
ergy spectrum of the incident solar neutrinos. 

The total event rate in the water experiments, 
about 0.5 the standard model value (see Fig. ||), 
is determined by the same high-energy 8 B neutri- 
nos that are expected, on the basis of the com- 
bined standard model, to dominate the event rate 
in the chlorine experiment. I have shown else- 
where that solar physics changes the shape 
of the 8 B neutrino spectrum by less than 1 part 
in 10 5 . Therefore, we can calculate the rate in 
the chlorine experiment (threshold 0.8 MeV) that 
is produced by the 8 B neutrinos observed in the 
Kamiokande and SuperKamiokande experiments 
at an order of magnitude higher energy threshold. 

If no new physics changes the shape of the 
8 B neutrino energy spectrum, the chlorine rate 
from 8 B alone is 2.8 ± 0.1 SNU for the Super- 
Kamiokande normalization (3.2 ±0.4 SNU for the 
Kamiokande normalization), which exceeds the 
total observed chlorine rate of 2.56 ± 0.23 SNU. 

Comparing the rates of the SuperKamiokande 
and the chlorine experiments, one finds-assuming 
that the shape of the energy spectrum of 8 B 
z/ e 's is not changed by new physics-that the net 
contribution to the chlorine experiment from the 
pep, 7 Be, and CNO neutrino sources is negative: 
—0.2 ± 0.3 SNU. The contributions from the pep, 
7 Be, and CNO neutrinos would appear to be com- 
pletely missing; the standard model prediction for 
the combined contribution of pep, 7 Be, and CNO 
neutrinos is a relatively large 1.8 SNU (see Ta- 
ble |l|) . On the other hand, we know that the 7 Be 
neutrinos must be created in the sun since they 
are produced by electron capture on the same iso- 
tope ( 7 Be) which gives rise to the 8 B neutrinos by 
proton capture. 

Hans Bethe and I pointed out [^8| that this ap- 
parent incompatibility of the chlorine and water- 
Cerenkov experiments constitutes a "second" so- 
lar neutrino problem that is almost independent 
of the absolute rates predicted by solar models. 
The inference that is usually made from this com- 
parison is that the energy spectrum of 8 B neu- 
trinos is changed from the standard shape by 
physics not included in the simplest version of 
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the standard electroweak model. 

3.3. Gallium Experiments: No Room for 
7 Be Neutrinos 

The results of the gallium experiments, 
GALLEX and SAGE, constitute the third solar 
neutrino problem. The average observed rate 
in these two experiments is 73 ± 5 SNU, which 
is accounted for in the standard model by the 
theoretical rate of 72.4 SNU that is calculated 
to come from the basic p-p and pep neutrinos 
(with only a 1% uncertainty in the standard solar 
model p-p flux). The 8 B neutrinos, which are ob- 
served above 6.5 MeV in the Kamiokande exper- 
iment, must also contribute to the gallium event 
rate. Using the standard shape for the spectrum 
of 8 B neutrinos and normalizing to the rate ob- 
served in Kamiokande, 8 B contributes another 6 
SNU. (The contribution predicted by the stan- 
dard model is 12 SNU, see Table @.) Given the 
measured rates in the gallium experiments, there 
is no room for the additional 34 ± 3 SNU that 
is expected from 7 Be neutrinos on the basis of 
standard solar models (see Table |l|) . 

The seeming exclusion of everything but p- 
p neutrinos in the gallium experiments is the 
"third" solar neutrino problem. This prob- 
lem is essentially independent of the previously- 
discussed solar neutrino problems, since it de- 
pends strongly upon the p-p neutrinos that are 
not observed in the other experiments and whose 
theoretical flux can be calculated accurately. 

The missing 7 Be neutrinos cannot be explained 
away by a change in solar physics. The 8 B neu- 
trinos that are observed in the Kamiokande ex- 
periment are produced in competition with the 
missing 7 Bc neutrinos; the competition is between 
electron capture on 7 Be versus proton capture on 
7 Be. Solar model explanations that reduce the 
predicted 7 Be flux generically reduce much more 
(too much) the predictions for the observed 8 B 
flux. 

The flux of 7 Be neutrinos, 0( 7 Be), is indepen- 
dent of measurement uncertainties in the cross 
section for the nuclear reaction 7 Be(p, 7) 8 B; the 
cross section for this proton-capture reaction is 
the most uncertain quantity that enters in an im- 
portant way in the solar model calculations. The 



flux of 7 Be neutrinos depends upon the proton- 
capture reaction only through the ratio 



( 7 Be) 



R(e) 



R(e) + R{p) ' 



(2) 



where R(e) is the rate of electron capture by 7 Be 
nuclei and R(p) is the rate of proton capture by 
7 Be. With standard parameters, solar models 
yield R(p) ~ 10~ 3 R(e). Therefore, one would 
have to increase the value of the 7 Be(p, 7) 8 B cross 
section by more than two orders of magnitude 
over the current best-estimate (which has an esti- 
mated experimental uncertainty of ~ 10%) in or- 
der to affect significantly the calculated 7 Be solar 
neutrino flux. The required change in the nuclear 
physics cross section would also increase the pre- 
dicted neutrino event rate by more than 100 in the 
Kamiokande experiment, making that prediction 
completely inconsistent with what is observed. 

I conclude that either: 1) at least three of the 
five operating solar neutrino experiments (the two 
gallium experiments plus either chlorine or the 
two water Cerenkov experiments, Kamiokande 
and SuperKamiokande) have yielded misleading 
results, or 2) physics beyond the standard elec- 
troweak model is required to change the energy 
spectrum of v e after the neutrinos are produced 
in the center of the sun. 

4. Uncertainties in the Flux Calculations 

I will now discuss uncertainties in the solar 
model flux calculations. 

Table ^| summarizes the uncertainties in the 
most important solar neutrino fluxes and in the 
CI and Ga event rates due to different nuclear fu- 
sion reactions (the first four entries), the heavy 
element to hydrogen mass ratio (Z/X), the ra- 
diative opacity, the solar luminosity, the assumed 
solar age, and the helium and heavy element dif- 
fusion coefficients. The 14 N +p reaction causes a 
0.2% uncertainty in the predicted pp flux and a 
0.1 SNU uncertainty in the CI (Ga) event rates. 

The predicted event rates for the chlorine and 
gallium experiments use recent improved calcula- 
tions of neutrino absorption cross sections. | p7| , p^| 
The uncertainty in the prediction for the gallium 
rate is dominated by uncertainties in the neutrino 
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Table 2 

Average uncertainties in neutrino fluxes and event rates due to different input data. The flux uncertainties 
are expressed in fractions of the total flux and the event rate uncertainties are expressed in SNU. The 
7 Be electron capture rate causes an uncertainty of ±2% |^9| that affects only the 7 Be neutrino flux. The 
average fractional uncertainties for individual parameters are shown. 



< Fractional 


PP 


3 He 3 He 


3 He 4 He 


7 Be +p 


Z/X 


opac 


lum 


age 


diffuse 


uncertainty > 


0.017 


0.060 


0.094 


0.106 


0.033 


see text 


0.004 


0.004 


0.15 


Flux 




















PP 


0.002 


0.002 


0.005 


0.000 


0.002 


0.003 


0.003 


0.0 


0.003 


7 Be 


0.0155 


0.023 


0.080 


0.000 


0.019 


0.028 


0.014 


0.003 


0.018 


8 B 


0.040 


0.021 


0.075 


0.105 


0.042 


0.052 


0.028 


0.006 


0.040 


SNUs 




















CI 


0.3 


0.2 


0.5 


0.6 


0.3 


0.4 


0.2 


0.04 


0.3 


Ga 


1.3 


0.9 


3.3 


1.3 


1.6 


1.8 


1.3 


0.20 


1.5 



absorption cross sections, +6.7 SNU (7% of the 
predicted rate) and —3.8 SNU (3% of the pre- 
dicted rate). The uncertainties in the chlorine ab- 
sorption cross sections cause an error, ±0.2 SNU 
(3% of the predicted rate), that is relatively small 
compared to other uncertainties in predicting the 
rate for this experiment. For non-standard neu- 
trino energy spectra that result from new neu- 
trino physics, the uncertainties in the predictions 
for currently favored solutions (which reduce the 
contributions from the least well-determined 8 B 
neutrinos) will in general be less than the val- 
ues quoted here for standard spectra and must 
be calculated using the appropriate cross section 
uncertainty for each neutrino energy. JlTjJl^] 

The nuclear fusion uncertainties in Table ^ 
were taken from Adelberger et al., [Q the neu- 
trino cross section uncertainties from Bahcall 
(1997)0 and Bahcall et al. (1996), @ the 
heavy element uncertainty was taken from he- 
lioseismological measurements, (3(]] the luminos- 
ity and age uncertainties were adopted from 
BP95, the lcr fractional uncertainty in the 
diffusion rate was taken to be 15%, J31J which 
is supported by helioseismological evidence, |23|] 
and the opacity uncertainty was determined by 
comparing the results of fluxes computed using 
the older Los Alamos opacities with fluxes com- 
puted using the modern Livermore opacities. |24|] 



To include the effects of asymmetric errors, the 
now publicly-available code for calculating rates 
and uncertainties (see discussion in previous sec- 
tion) was run with different input uncertainties 
and the results averaged. The software contains a 
description of how each of the uncertainties listed 
in Table || were determined and used. 

The low energy cross section of the 7 Be + p re- 
action is the most important quantity that must 
be determined more accurately in order to de- 
crease the error in the predicted event rates in 
solar neutrino experiments. The 8 B neutrino flux 
that is measured by the Kamiokande, || Super- 
Kamiokande, || and SNO J32j experiments is, 
in all standard solar model calculations, directly 
proportional to the 7 Be+p cross section. If the la 
uncertainty in this cross section can be reduced 
by a factor of two to 5%, then it will no longer be 
the limiting uncertainty in predicting the crucial 
8 B neutrino flux (cf. Table |). 

5. How Large an Uncertainty Does Helio- 
seismology Suggest? 

Could the solar model calculations be wrong by 
enough to explain the discrepancies between pre- 
dictions and measurements for solar neutrino ex- 
periments? Helioseismology, which confirms pre- 
dictions of the standard solar model to high pre- 
cision, suggests that the answer is probably "No." 
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Figure ^ shows the fractional differences be- 
tween the most accurate available sound speeds 
measured by helioseismology f33| and sound 
speeds calculated with our best solar model (with 
no free parameters). The horizontal line cor- 
responds to the hypothetical case in which the 
model predictions exactly match the observed 
values. The rms fractional difference between 
the calculated and the measured sound speeds 
is 1.1 x 1CP 3 for the entire region over which 
the sound speeds are measured, 0.05i? Q < R < 
0.95i? Q . In the solar core, 0.05i? Q < R < 0.25i? Q 
(in which about 95% of the solar energy and neu- 
trino flux is produced in a standard model), the 
rms fractional difference between measured and 
calculated sound speeds is 0.7 x 10~ 3 . 



o.i 

0.08 

0.06 

S 0.04 
m 

> 0.02 
c 

m 

| -0.02 
o 

5 -0.04 
-0.06 
-0.08 
-0.1 





I i i i I 




i i i i i i i i 






Bahcall 


-Pinsonneault 98 






L0WL1 - 


- BISON Measurement - 


- <- 3a 


v discrepancy 






■ «- la 


v discrepancy 








i . , . i 







0.2 0.4 0.6 0.8 1 

R/R s 

Figure 6. Predicted versus Measured Sound 
Speeds. This figure shows the excellent agree- 
ment between the calculated (solar model BP98, 
Model) and the measured (Sun) sound speeds, a 
fractional difference of 0.001 rms for all speeds 
measured between 0.05i?o and 0.95i?Q. The 
vertical scale is chosen so as to emphasize that 
the fractional error is much smaller than generic 
changes in the model, 0.04 to 0.09, that might 
significantly affect the solar neutrino predictions. 



Helioseismological measurements also deter- 
mine two other parameters that help character- 



ize the outer part of the sun (far from the in- 
ner region in which neutrinos are produced): the 
depth of the solar convective zone (CZ), the re- 
gion in the outer part of the sun that is fully 
convective, and the present-day surface abun- 
dance by mass of helium (Y" SU rf)- The mea- 
sured values, R C z = (0.713 ± 0.001)i? Q , @ and 
F surf = 0.249 ± 0.003, [|0) are in satisfactory 
agreement with the values predicted by the so- 
lar model BP98, namely, Rcz = 0.714i? Q , and 
Y sur f = 0.243. However, we shall see below that 
precision measurements of the sound speed near 
the transition between the radiative interior (in 
which energy is transported by radiation) and the 
outer convective zone (in which energy is trans- 
ported by convection) reveal small discrepancies 
between the model predictions and the observa- 
tions in this region. 

If solar physics were responsible for the solar 
neutrino problems, how large would one expect 
the discrepancies to be between solar model pre- 
dictions and helioseismological observations? The 
characteristic size of the discrepancies can be es- 
timated using the results of the neutrino exper- 
iments and scaling laws for neutrino fluxes and 
sound speeds. 

All recently published solar models predict es- 
sentially the same fluxes from the fundamental 
pp and pep reactions (amounting to 72.4 SNU 
in gallium experiments, cf. Table 1), which are 
closely related to the solar luminosity. Compar- 
ing the measured gallium rates and the standard 
predicted rate for the gallium experiments, the 
7 Be flux must be reduced by a factor N if the dis- 
agreement is not to exceed n standard deviations, 
where N and n satisfy 72.4 + (34.4) /iV = 72.2 + 
na. For a la (3a) disagreement, N = 6.1(2.05). 
Sound speeds scale like the square root of the lo- 
cal temperature divided by the mean molecular 
weight and the 7 Be neutrino flux scales approxi- 
mately as the 10th power of the temperature. |35| 
Assuming that the temperature changes are dom- 
inant, agreement to within la would require frac- 
tional changes of order 0.09 in sound speeds (3a 
could be reached with 0.04 changes), if all model 
changes were in the temperature^. This argu- 



1 I have used in this calculation the GALLEX and SAGE 
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ment is conservative because it ignores the 8 B and 
CNO neutrinos which contribute to the observed 
counting rate (cf. Table 1) and which, if included, 
would require an even larger reduction of the Be 
flux. 

I have chosen the vertical scale in Fig. ^ to 
be appropriate for fractional differences between 
measured and predicted sound speeds that are of 
order 0.04 to 0.09 and that might therefore af- 
fect solar neutrino calculations. Fig. ^ shows that 
the characteristic agreement between solar model 
predictions and helioseismological measurements 
is more than a factor of 40 better than would be 
expected if there were a solar model explanation 
of the solar neutrino problems. 

6. Fits Without Solar Models 

Suppose (following the precepts of Hata et 
al., [3q] Parke, J3?J and Heeger and Robert- 
son [38|) we now ignore everything we have 
learned about solar models over the last 35 years 
and allow the important pp, 7 Be, and 8 B fluxes 
to take on any non- negative values. What is the 
best fit that one can obtain to the solar neutrino 
measurements assuming only that the luminosity 
of the sun is supplied by nuclear fusion reactions 
among light elements (the so-called 'luminosity 
constraint')? Q 

The answer is that the fits are bad, even if 
we completely ignore what we know about the 
sun. I quote the results from Bahcall, Krastev 
and Smirnov (1998). 

If the CNO neutrino fluxes are set equal to zero, 
there are no acceptable solutions at the 99% C. 
L. (~ 3a result). The best-fit is worse if the CNO 
fluxes are not set equal to zero. All so-called 'so- 
lutions' of the solar neutrino problems in which 
the astrophysical model is changed arbitrarily (ig- 
noring helioseismology and other constraints) are 
inconsistent with the observations at much more 

measured rates reported by Kirsten and Gavrin at Neu- 
trino 98. The experimental rates used in BP98 were not 
as precise and therefore resulted in slightly less stringent 
constraints than those imposed here. In BP98, we found 
that agreement to within la with the then available ex- 
perimental numbers would require fractional changes of 
order 0.08 in sound speeds (3<r could be reached with 0.03 
changes.) 



Table 3 

Neutrino Oscillation Solutions. 



Solution 


Am 2 


sin 2 29 


SMA 


5 x 10~ 6 eV 2 


5 x 10" 3 


LMA 


2 x 10" 5 eV 2 


0.8 


LOW 


8 x 10~ 8 eV 2 


0.96 


VAC 


8 x 10~ n eV 2 


0.7 



than a 3a level of significance. No fiddling of 
the physical conditions in the model can yield the 
minimum value, quoted above, that was found by 
varying the fluxes independently and arbitrarily. 

Figure ^ shows, in the lower left-hand corner, 
the best-fit solution and the lcr -3a contours. The 
la and 3c limits were obtained by requiring that 
X 2 = Xmin + ^X 2 > where for la 8\ 2 — 1 and for 
3a 5x 2 = 9. All of the standard model solutions 
lie far from the best-fit solution and even lie far 
from the 3a contour. 

Since standard model descriptions do not fit the 
solar neutrino data, we will now consider models 
in which neutrino oscillations change the shape of 
the neutrino energy spectra. 

7. Neutrino Oscillations 

The experimental results from all five of 
the operating solar neutrino experiments (chlo- 
rine, Kamiokande, SAGE, GALLEX, and Super- 
Kamiokande) can be fit well by descriptions in- 
volving neutrino oscillations, either vacuum os- 
cillations (as originally suggested by Gribov and 
Pontecorvo |h|)or resonant matter oscillations 
(as originally discussed by Mikheyev, Smirnov, 
and Wolfenstein (MSW) jS)). 

Table || summarizes the four best-fit solu- 
tions that arc found in the two-neutrino ap- 
proximation. pl],p6l Only the SMA and vac- 
uum oscillation solutions fit well the recoil elec- 
tron energy spectrum measured in the Super- 
Kamiokande experiment-if the standard value 
for the hep production reaction cross section 
( 3 He + p -> 4 He + e+ + v e ) is used. Jll[] 
However, for over a decade I have not given an es- 
timated uncertainty for this cross section. B The 
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transition matrix element is essentially forbidden 
and the actual quoted value for the production 
cross section depends upon a delicate cancellation 
between two comparably sized terms that arise 
from very different and hard to evaluate nuclear 
physics. I do not see anyway at present to deter- 
mine from experiment or from first principles the- 
oretical calculations a relevant, robust upper limit 
to the hep production cross section (and therefore 
the hep solar neutrino flux). 

The possible role of hep neutrinos in solar neu- 
trino experiments is discussed extensively in Bah- 
call and Krastev (1998) |2(| The most important 
unsolved problem in theoretical nuclear physics 
related to solar neutrinos is the range of values 
allowed by fundamental physics for the hep pro- 
duction cross section. 

8. Discussion 

When the chlorine solar neutrino experiment 
was first proposed, the only stated motiva- 
tion was "...to see into the interior of a star and 
thus verify directly the hypothesis of nuclear en- 
ergy generation in stars." This goal has now been 
achieved, 

The focus has shifted to using solar neutrino 
experiments as a tool for learning more about the 
fundamental characteristics of neutrinos as parti- 
cles. Experimental effort is now concentrated on 
answering the question: What are the probabil- 
ities for transforming a solar v e of a definite en- 
ergy into the other possible neutrino states? Once 
this question is answered, we can calculate what 
happens to 24 's that are created in the interior 
of the sun. Armed with this information from 
weak interaction physics, we can return again to 
the original motivation of using neutrinos to make 
detailed, quantitative tests of nuclear fusion rates 
in the solar interior. Measurements of the fla- 
vor content of the dominant low energy neutrino 
sources, p-p and 7 Be neutrinos, will be crucial 
in this endeavor and will require another genera- 
tion of superb solar neutrino experiments (see the 
comments in Section ^) . 

Three decades of refining the input data and 
the solar model calculations has led to a predicted 
standard model event rate for the chlorine exper- 



iment, 7.7 SNU, which is very close to 7.5 SNU, 
the best-estimate value obtained in 1968. § The 
situation regarding solar neutrinos is, however, 
completely different now, thirty years later. Four 
experiments have confirmed the original chlorine 
detection of solar neutrinos. Helioseismological 
measurements are in excellent agreement with 
the standard solar model predictions and very 
strongly disfavor (by a factor of 40 or more) hypo- 
thetical deviations from the standard model that 
are require to fits the neutrino data (cf. Fig. ||). 
Just in the last two years, improvements in the 
helioseismological measurements have resulted in 
a five-fold improvement in the agreement between 
the calculated standard solar model sound speeds 
and the measured solar velocities (cf. Figure 2 of 
the Neutrino 96 talk |@ with Figure | of this 
talk). 

9. What next? 

More than 98% of the calculated standard 
model solar neutrino flux lies below 1 MeV. The 
rare 8 B neutrino flux is the only solar neutrino 
source for which measurements of the energy have 
been made, but 8 B neutrinos constitute a fraction 
of less than 10 -4 of the total solar neutrino flux. 

The next goal of solar neutrino astronomy is to 
measure neutrino fluxes below 1 MeV. We should 
begin today preparing for experiments that will 
measure the 7 Be neutrinos (energy of 0.86 MeV) 
and the fundamental p-p neutrinos (< 0.43 MeV). 
Indeed, we have heard at this workshop some 
marvelously exciting descriptions of how such low 
energy experiments could be carried out. The 
BOREXINO observatory, which can detect v — e 
scattering, is the only approved solar neutrino ex- 
periment which can measure energies less than 1 
MeV. 

The p-p neutrinos are overwhelmingly the most 
abundant source of solar neutrinos, carrying 
about 91% of the total flux according to the stan- 
dard solar model. The 7 Be neutrinos constitute 
about 7% of the total standard model flux. 

If we want to test and to understand neu- 
trino oscillations with high precision using solar 
neutrino sources, then we have to measure the 
neutrino-electron scattering rate with 7 Be neutri- 
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nos, as will be done with the BOREXINO experi- 
ment, and also the CC (neutrino-absorption) rate 
with 7 Be neutrinos (no approved experiment). 
With a neutrino line as provided by 7 Be electron- 
capture in the sun, unique and unambiguous tests 
of neutrino oscillation models can be carried out 
if one knows both the charged-current and the 
neutral current reaction rates [Q . 

I believe we have calculated the flux of p-p neu- 
trinos produced in the sun to an accuracy of ±1%. 
Unfortunately, we do not yet have a direct mea- 
surement of this flux. The gallium experiments 
only tell us the rate of capture of all neutrinos 
with energies above 0.23 MeV. 

The most urgent need for solar neutrino re- 
search is to develop a practical experiment to 
measure directly the p-p neutrino flux and the 
energy spectrum of electrons produced by target 
interactions with p-p neutrinos. Such an experi- 
ment can be used to test the precise and funda- 
mental standard solar model prediction of the p- 
p neutrino flux. Moreover, the currently favored 
neutrino oscillation solutions all predict a strong 
influence of oscillations on the low-energy flux of 
v e . 

Figure |?] shows the calculated neutrino survival 
probability as a function of energy for three global 
best-fit MSW oscillation solutions. You can see 
directly from this figure why we have to have ac- 
curate measurements for the p-p and 7 Be neutri- 
nos: the currently favored solutions exhibit their 
most characteristic and strongly energy depen- 
dent features below 1 MeV. In all of these solu- 
tions, the survival probability shows a dramatic 
increase with energy below 1 MeV, whereas in 
the region above 5 MeV (accessible to Super- 
Kamiokande and to SNO) the energy dependence 
of the survival probability is at best modest. 

The p-p neutrinos are the gold ring of solar neu- 
trino astronomy. Their measurement will con- 
stitute a simultaneous and critical test of stellar 
evolution theory and of neutrino oscillation solu- 
tions. 

The most exciting result of this workshop for 
me has been the possibility discussed here of a 
synergistic experiment involving a huge (mega- 
ton?) nucleon decay detector in which an in- 
ner region is reserved for solar neutrino experi- 
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Figure 7. Survival probabilities for MSW solu- 
tions. The figure presents the yearly-averaged 
survival probabilities for an electron neutrino that 
is created in the sun to remain an electron neu- 
trino upon arrival at the SuperKamiokande de- 
tector. There are only slight differences be- 
tween the computed regeneration probabilities for 
the detectors located at the positions of Super- 
Kamiokande, SNO and the Gran Sasso Under- 
ground Laboratory. The full line refers to the av- 
erage survival probabilities computed taking into 
account regeneration in the earth and the dotted 
line refers to calculations for the daytime that 
do not include regeneration. The dashed line 
includes regeneration at night. This is Fig. 9 
of the 1998 analysis by Bahcall, Krastev, and 
Smirnov pj]. 



ments(see, for example, the talks by Jung, Naka- 
hata, and Ypsilantis at this workshop) . The most 
straightforward solar neutrino experiments that 
could be carried out with this detector would be 
precision measurements of the temporal depen- 
dences of the relatively high-energy 8 B neutri- 
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nos. One could measure with such a detector 
the zenith-angle dependence of the solar neutrino- 
event rate (the generalization of the day-night dif- 
ference) and the seasonal dependence (generaliza- 
tion of the winter-summer difference) . The design 
could relax somewhat the precise requirements for 
energy calibration and for energy resolution used 
for the SuperKamiokande and SNO experiments 
and concentrate instead on limiting the system- 
atic uncertainties in the detector that could con- 
tribute to the error budget in the day-night or 
seasonal dependences. After three years of very 
careful measurements, the SuperKamiokande ex- 
periment has, as we have heard at this conference, 
about a 2a result for the day-night difference. 
They do not yet have the statistics to report a 
meaningful measurement of the full zenith-angle 
dependence or the seasonal dependence. The pre- 
dicted temporal effects are small, generally of or- 
der a percent, with the currently favored neutrino 
oscillation solutions. 

Nature has provided us with many different 
baselines and with many different matter column 
densities with which to do Very Long Baseline 
(VLB) studies of neutrino oscillations. The earth- 
sun distance varies continuously during the year 
between 1.496(1.0 ± 0.017)10 13 cm and the col- 
umn density through the earth to a terrestrial 
detector varies from gm cm~ 2 during the day 
to more than 10 9 gm cm -2 at night. 

A solar neutrino detector ten or more times the 
volume of the current SuperKamiokande experi- 
ment, as discussed in concept at this workshop, 
could measure precisely the results of many differ- 
ent VLB neutrino oscillation experiments. This 
would be a fantastic 'Smoking Gun' detector. I 
have a hard time sitting down when imagining 
such an exciting possibility. 
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